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Abstract 

We analyze eigenvalues fluctuations of the Laplacian of various networks under 
the random matrix theory framework. Analyses of random networks, scale- free 
networks and small-world networks show that nearest neighbor spacing distri- 
bution of the Laplacian of these networks follow Gaussian orthogonal ensemble 
statistics of random matrix theory. Furthermore, we study nearest neighbor spac- 
ing distribution as a function of the random connections and find that transition to 
the Gaussian orthogonal ensemble statistics occurs at the small-world transition. 
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1 Introduction 

In order to understand complex world around us network theory has been getting fast 
recognition. The main concept here is to define complex systems in terms of networks 
of many interacting units. Few examples of such systems are interacting molecules in 
living cell, nerve cells in brain, computers in Internet communication, social networks of 
interacting people, airport networks with flight connections, etc [H [21 El S] . To under- 
stand these networks simple models are introduced; these model networks are based on 
some simple principles, and capture essential features of real systems. Mathematically 
networks are investigated under the framework of graph theory. In the graph theoretical 
terminology, units are called nodes and interactions are called edges [5]. 

Random graph model of Erdos and Renyi (ER) assumes that interaction between 
the nodes are random [6]. Recently, with the availability of large maps of real world 
networks, it has been observed that the random graph model is not appropriate for 
studying real world networks. Hence many new models have been introduced. Barabasi- 
Albert's scale-free (SF) model [7j and Watts-Strogatz's small-world (SW) model [8] 
are the most recognized ones, and have contributed immensely in understanding the 



evolution and behavior of the real systems having network structures. SF model captures 
the degree distribution behavior of real world networks and SW model deals with the 
clustering coefficient and diameter. Following these two new models came an outbreak 
in the field of networks. These works have focused on the following aspects: (1) direct 
studies of the real-world networks and measuring their various structural properties such 
as degree distribution, diameter, clustering coefficient, etc., (2) proposing new random 
graph models motivated by these studies, and (3) computer simulations of the new 
models and measuring their properties [2]. 

Barabasi et.al. investigations of the various real world systems show that they are 
scale-free, which means that the degree distribution p{k), fraction of nodes that have k 
number of connections with other nodes, decays as power law, i.e. p{k) oc k~'^, where 
7 depends on the topology of the networks. The scale-free nature of networks suggests 
that there exist few nodes with very high degree. Some other analysis, by Newman and 
others, of real-world networks show that complex networks have community or module 
structures [9l [10]. According to these studies, there exists few nodes with very high 
betweenness which are responsible to connect the different communities. This direction 
of looking at the networks focuses on the importance of nodes based on its position in the 
network. On the other hand, ER and SW models emphasize on the random connections 
in the networks; in ER model any two nodes are connected with probability p. One 
of the most interesting characteristics of ER model was sudden emergence of various 
global properties; most important one being emergence of giant cluster. For a. p > pc, 
while number of nodes in the graph remain constant, giant cluster emerges through a 
phase transition. Further, Watts-Strogatz model shows the small-world transition with 
fine tuning of the number of random connections. 

Apart from the above mentioned investigations which focus on direct measurements 
of structural properties of networks, there exists a vast literature demonstrating that 
properties of networks or graphs could be well characterized by the spectrum of associ- 
ated adjacency (A) and Laplacian (L) matrix [ll]. For an unweighted graph, adjacency 
matrix is defined in following way : Aij = 1, if i and j nodes are connected and zero 
otherwise. Laplacian of graph has been defined in the various ways (depending upon 
the normalization) in the literature. We follow definition used in [.12j; 



where di denotes the degree of node i. For undirected networks, adjacency and Laplacian 
both are symmetric matrices and consequently have real eigenvalues [13]. Eigenvalues 
of graph are called graph spectra and they give informations about some basic topolog- 
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ical properties of the underlying network. During last twenty years several important 
applications of the spectral graph theory in physics and chemistry problems have been 
discovered [TTl [T4] . For example liquid flowing through a system of communicating pipes 
are described by a system of linear differential equations. The corresponding matrix ap- 
pears to be the Laplacian of the underlying graph. Speed of convergence of the liquid 
flowing process towards an equilibrium state is measured by the second smallest eigen- 
value of L [13]. Second smallest eigenvalue of L is also called the algebraic connectivity 
of a graph and is used to understand behavior of dynamical processes on the underlying 
networks [15]. Particularly, Laplacian spectra of networks have been investigated enor- 
mously to understand synchronization of coupled dynamics on networks [151 [El flTj - 
for example recently extremal eigenvalues of the Laplacian have been shown to have 
high influence on the synchronizability of the network [18]. Similarly, multiplicities of 
eigenvalues, particularly at and 1, have direct relations with the properties of graphs 

[iniEo]. 

Following our recent works, where we used random matrix theory (RMT) to study 
spectral properties of adjacency matrix of various networks [2T], in this paper we investi- 
gate spectral properties of Laplacians of networks under RMT framework. Particularly 
we study nearest neighbor spacing distribution (NNSD) of Laplacian matrix of various 
model networks, namely scale-free, small world and random networks. We find that 
inspite of spectral densities of different model networks are different, their eigenvalue 
fluctuations are same and follow Gaussian orthogonal ensemble (GOE) distribution of 
RMT. We attribute this universality to the presence of the similar amount of random- 
ness in all these networks, and show that randomness in the network connections can be 
quantified by the Brody parameter coming from RMT. Furthermore, there exists one to 
one correlation between the diameter of the network and the eigenvalues fluctuations of 
the Laplacian matrix. By changing number of connections in the network we get tran- 
sition to the GOE distribution. As Erdos and Renyi observed that with the fine tuning 
of network parameter all nodes get connected with a sudden transition; under the RMT 
framework our analysis suggests transition to some kind of spreading of randomness 
over the whole network. 

Note that in this paper we consider normalized Laplacians, though for the RMT 
analysis form of Laplacians does not matter, because unfolding, a method to separate out 
system dependent part from the eigenvalues to study their universal behavior, removes 
scaling effects caused by the normalization factor. We consider normalized Laplacians 
because eigenvalues properties of this form of Laplacian are extensively investigated 
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[22| [231 1211 [251 126]. They have been found to be an excellent candidate as a concise 
fingerprint of internet-like graphs [2lj and graphs for the cancer cells [25] . Furthermore 
for some applications, such as graph partitioning, this normalized form is preferred over 
other forms of Laplacian [261 [20] - 

The paper is organized as follows: after this introductory section, in Sec. [2l we 
describe some techniques of RMT which we have used in our analysis. In Sec. [3l we 
analyze the NNSD of Laplacian for various networks, namely small- world, scale-free and 
random networks. In section [H we study the effect of random connections on the level 
statistics. In this section we show the one to one correlation between GOE transition 
and small-world behavior. Finally, in Sec. [3, we summarize and discuss about some 
possible future directions. 

2 Random matrix statistics 

RMT was proposed by Wigner to explain the statistical properties of nuclear spectra 
[27] . Later this theory had successfully been applied in the study of different complex 
systems such as disordered systems, quantum chaotic systems where RMT tells whether 
corresponding classical system is regular or chaotic or a mixture of both, spectra of large 
complex atoms, etc |28j. RMT is also shown to be of great interest in understanding 
the statistical structure of the empirical cross-correlation matrices appearing in the 
study of multivariate time series. The classical complex systems where RMT has been 
successfully applied are stock market (cross-correlation matrix is formed by using the 
time series of price fluctuations of different stock) [29]; brain (matrix is constructed by 
using EEG data at different locations) [30]; patterns of atmospheric variability (cross- 
correlations matrix is generated by using temporal variation of various atmospheric 
parameters) [3lj, etc. 

We study eigenvalues fluctuations of the Laplacian of the networks. The eigenval- 
ues fluctuations are generally obtained from the nearest neighbor spacing distribution 
(NNSD) of the eigenvalues. The NNSD follows two universal properties depending 
upon the underlying correlations among the eigenvalues. For correlated eigenvalues, the 
NNSD follows Wigner-Dyson formula of Gaussian orthogonal ensemble (GOE) statis- 
tics of RMT; whereas, the NNSD follows Poisson statistics of RMT for uncorrelated 
eigenvalues. 

Here we briefly describe some aspects of RMT which we have used in our network 
analysis. We denote the eigenvalues of network Laplacian by Aj, i = 1, . . . , N, where N 
is the size of the network and Aj+i > Aj, Vi. In order to get universal properties of the 
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fluctuations of tlie eigenvalues, it is customary in RMT to unfold the eigenvalues by a 
transformation Aj = N{Xi), where N is the averaged integrated eigenvalue density [27] . 
Since we do not have any analytical form for A^, we numerically unfold the spectrum 
by polynomial curve fitting. After the unfolding, the average spacing will be unity, 
independent of the system. Using the unfolded spectra, we calculate the spacing as 
Si = Aj+i — Aj. The NNSD P{s) is defined as the probability distribution of these Sj's. 
In case of Poisson statistics, P{s) = exp(— s); whereas for GOE 

P(s) = -sexpi^-—j. (2) 

For the intermediate cases, the spacing distribution is described by Brody distribution 

Pp{s) = As^ exp , (3) 

where 

13+1 

(4) 

This is a semiempirical formula characterized by the parameter j3. As (3 goes from to 
1, the Brody distribution smoothly changes from Poisson to GOE. We fit the spacing 
distributions of different networks by Brody distribution Pj3{s). This fitting gives an 
estimation of (3, and consequently identifies whether the spacing distribution of a given 
network is Poisson or GOE or intermediate of these two. 

3 Laplacian matrix spectrum of complex networks 

Following we present the ensemble averaged spectral density and spacing distribution 
of random, scale-free and small-world networks. 

3.1 Random network 

First we consider random network generated by using Erdos and Renyi algorithm. We 
take = 2000 nodes and make random connections between pairs of nodes with proba- 
bility p = 0.01. This method yields a connected network with average degree px A^ = 20. 
Note that for very small value of p one gets several unconnected components. Here, the 
choice of p is such that it should be high enough {p > pc) to give large connected com- 
ponent typically spanning all nodes [33]; and since most of real world networks are very 
sparse [2], p is small enough to have a sparse network. Laplacian is constructed using 
Eq. ([1]). Because of the normalized form, the eigenvalues would always be within and 
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Figure 1: (Color online) (a) and (b) show the spectral density (p(A)) and corresponding 
spacing distribution (-P(s)) for random network. The histograms are numerical results 
and the solid lines represent fitted Brody distribution. Figures are plotted for average 
over 10 random realizations of the networks. All networks have = 2000 nodes and 
an average degree k = 20 per node. 

2. Figure [T](a) plots the spectral density of the Laplacian. Note that the distribution is 
averaged over 10 realizations of the network. The spectral density follows Wigner-Dyson 
semicircular distribution [3l] . To get the spacing behaviors, first the eigenvalues are un- 
folded by using the technique described in Sec. [2l This method yields the eigenvalues 
with constant spectral density on the average. These unfolded eigenvalues are used to 
calculate NNSD. The same procedure is repeated for an ensemble of the networks gen- 
erated for different random realizations. Figure [11(b) plots ensemble average of NNSD. 
By fitting this ensemble averaged NNSD with the Brody formula given in Eq. ([3]) we get 
an estimation of the Brody parameter (3 = 0.9786 ~ 1. This value of Brody parameter 
clearly indicates the GOE behavior of the NNSD [Eq. (^]. 

3.2 Scale-free network 

Figure [2] plots the spectral density and the nearest neighbor spacing distribution of the 
Laplacian of scale-free network. Scale-free network is generated by using the model of 
Barabasi et al. [7]. Starting with a small number, mo of the nodes, a new node with 
m < niQ connections is added at each time step. This new node connects with node i 
with the probability 7r(fcj) oc ki (preferential attachment), where ki is the degree of the 
node i. After r time steps the model leads to a network with A^ = T + mQ nodes and mr 
connections. This model leads to a scale-free network, i.e., the probability P{k) that a 
node has degree k decays as a power law, P{k) ~ where A is a constant and for the 
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Figure 2: (Color online) (a) and (b) show the spectral density (p(A)) and corresponding 
spacing distribution (-P(s)) of scale-free network. The histograms are numerical results 
and the solid lines represent fitted Brody distribution. Figures are plotted for average 
over 10 random realizations of the networks. All networks have = 2000 nodes and 
an average degree k = 20 per node. 



type of probability law TT{k) that we have used A = 3. Other forms for the probability 
iT{k) are possible which give different values of A and we find results similar to the ones 
reported here. Density distribution of the network has a very uniform distribution for 
almost all eigenvalues except for few lowest and highest ones. There is a slight peak 
around one which corresponds to the peak at zero for the adjacency matrix (|21j). 
To calculate NNSD, we follow same procedure as described in the previous section. 
Fig. [2](b) shows that spacing distribution of this scale-free network follow GOE very 
closely {p = 0.9914 ~ 1). 

3.3 Small-world network 

Figure [3] is plotted for small- world network. Small- world networks are constructed using 
the following algorithm by Watts and Strogatz [8]. Starting with a one-dimension ring 
lattice of nodes in which every node is connected to its k nearest neighbors, we 
randomly rewire each connection of the lattice with the probability p such that self- 
loop and multiple connections are excluded. Thus p = gives a regular network and 
p = 1 gives a random network. The typical small world behavior is observed around 
p = 0.005. We take A^ = 2000 and average degree k = 40. Spectral density of this 
network is complicated with several peaks. One peak is at A = 1 which corresponds to 
the p{0) peak for the spectra of adjacency matrix. For different values of k the exact 
positions of other peaks may vary but overall form of spectral density remains similar. 



7 



Figure 3: (Color online) (a) and (b) show the spectral density (p(A)) and correspond- 
ing spacing distribution (-P(s)) of small-world network. The histograms are numerical 
results and the solid lines represent fitted Brody distribution. Figures are plotted for 
average over 10 random realizations of the networks. All networks have = 2000 nodes 
and an average degree k = 20 per node. 

Fig. Eb) shows that spacing distribution of small world network follows GOE very closely 
{p = 0.9584 ~ 1). 

Figs. [11(b), [2](b) and [3b) show that spacing distribution of all the three networks 
follow GOE very closely (/5 ^ 1). Since spectral density of random network is very close 
to the semicircular, and in RMT literature semicircular distributions are extensively 
studied and are shown to give GOE statistics of the spacings, we expected that NNSD of 
random networks would follow GOE as well. But for scale-free network and small- world 
network, density distributions have very different forms, still NNSD of these networks 
follow GOE statistics. This is a very interesting result. Following RMT these results 
imply that even though the spectral density of the scale-free and small-world network 
differ than the spectral density of the random network, but the correlations among the 
eigenvalues are strong enough to yield GOE statistics. 

4 Transition to GOE statistics and small-world be- 
havior 

All the above networks have some amount of random connections in them. In this 
section we discuss the networks which are regular and analyze the spectral properties 
as we go from regular network to the random one. Starting with a ring lattice of 
nodes A^ = 2000 and average degree k = 40, we then randomize some connections with 



8 



7 
6 

51- 
24- 
^3 

2 

1 





(a) 



. ' 1 ' 


' 1 ' . 
(b) : 




i 1 , ■ 





' 1 ' . 
(c) : 








Figure 4: Figures shows the transition from Poisson to the GOE statistics of the eigen- 
values of the Laplacian of the network as a function of random rewiring parameter p. 
We take regular lattice of size = 2000 and average degree k = 40. Links are rewires 
with probability p. Figures are plotted for average over 10 random realization of the 
networks. 



probability p, and observe the change in the spectral density and NNSD of the network 
as a function of p. For the original ring lattice, the underlying adjacency matrix would 
be a band matrix with entries one in the band, except diagonal elements which would 
be zero. Corresponding Laplacian would be a band matrix as well, with off-diagonal 
elements x k (from Eq. ([1])), where k being the degree of each node, and diagonal 

elements 1. As a result of rewiring, few randomly selected pairs of nodes get connected, 
yielding few nonzero entries outside to the band of the corresponding Laplacian, while 
making few band elements zero. For a completely random network, the corresponding 
Laplacian would have non zero random entries. 

Following we discuss the spectral behavior of Laplacian matrix as number of ran- 
dom connections in the underlying network is increased. Figures H] (a) and (d) plot 
density distribution and NNSD for the p = 5 x 10~^ case. At this value of p there is 
as few as one random rewiring in the regular lattice. These figures shows that spectral 
density of the lattice is complicated without having any known form; and its spacing 
distribution closely follows Poisson statistics (/5 ~ 0.1). We then increase the value of p 
and randomize a fraction p = 2x 10~^ of the edges. For this value of p, the spectral den- 
sity and the spacing distribution are plotted respectively in Figure |l](b) and (e). These 
figures reveal that for this small value of p, density distribution of the network does 
not show any noticeable change, whereas spacing distribution shows different property 
(corresponding to the P ~ 0.67). As p is further increased to 5 x 10""^ the spectral den- 
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sity hardly shows any change but interestingly spacing distribution shows significantly 
different property than the Poisson statistics. Now the spacing distribution looks more 
like the intermediate between Poisson and GOE. By fitting the spacing distribution cor- 
responding to this p value with the Brody formula, we estimate the Brody parameter 
/3 as /3 ~ 0.8. As marked change in the Brody parameter for very small change in the 
random connections in the network (corresponding to p = 5 x 10~^, Fig. H] (d)), we had 
expected that the Brody parameter will reach asymptotically to unity as we increase p. 
But the surprising finding is that the onset of transition from Poisson to GOE occurs 
at very small value of p (Figs. H] (e) and (f)), and is related with the SW properties of 
the networks. We calculate the network diameters and clustering coefficients for these 
values of p, which are listed below: 

Table 1: Diameter, clustering coefficient, and Brody parameter for the various values of 
rewiring probability p. 





p 


m/m 


Cip)/C{0) 


1-/3 


(a) 


5 X 10^5 


0.83925 


0.99980 


0.89463 


(b) 


2 X 10"^ 


0.53680 


0.99925 


0.33368 


(c) 


5 X 10"^ 


0.39776 


0.99857 


0.20131 



Table [H shows that for (b) and (c), the average length of the network is as small as 
the corresponding random network {Lrandom ~ ln{N)/\n{k), and clustering coefficient 
is as high as the regular lattice {Creguiar ~ 3/4), which are the properties of small- world 
networks. Values of (3 reveals that Poisson to GOE (/3 ~ ^ /3 ~ 1) transition and 
small-world transition takes place for the similar value of p. 

5 Conclusions 

In summary, we study eigenvalues spacing distribution of Laplacian of the model net- 
works studied extensively in the literature. They follow universal GOE statistics. From 
RMT analogy it tells that there exists correlations between the eigenvalues of the net- 
work arising from certain symmetries in the interactions. We attribute this universality 
to the similar amount of randomness in the model networks, which may arise natu- 
rally, in order to capture the real world network properties. We study the effect of 
the randomness in the network connections on the eigenvalues fluctuations of network 
Laplacian, and use Brody parameter to quantify this randomness. These studies reveal 
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that there is a direct relation between the random connections and the Brody param- 
eter. For the regular network (with the average degree greater then some value) we 
get Poisson distribution, as we make random rewiring of the connections the spacing 
distribution start deviating from the Poisson statistics, first it shows intermediate of 
Poisson and GOE statistics and for sufficiently large number of random connections it 
show almost GOE statistics {(3 ~ 1). According to the interpretation of RMT, at this 
value of p, eigenvalues are as much correlated as for the completely random networks. 
Furthermore we find that this transition to GOE statistics happens at the onset of the 
small world behavior. 

Universal RMT results shown by the networks suggest that they can be modeled 
as a random matrix chosen from GOE of RMT. Eigenvalues fluctuations following GOE 
statistics argue that there exists some kind of spreading over the randomness in the 
whole network, which may be essential for robustness of the system. According to 
many recent studies, randomness in the connections is one of the most important and 
desirable ingredients for the proper functionality or the efficient performance of the 
system having underlying network structure. For instance, information processing in 
the brain is considered to be highly influenced by random connections among different 
modular structure [35] . Our analysis suggests that randomness in the complex networks 
can be studied under the RMT framework. Furthermore, Laplacian spectra have been 
investigated to understand various dynamical processes on networks |16l [171 [IE] as well, 
hence the RMT analysis of network Laplacians could also be important in view of these 
demonstrations of the relations between spectral properties and dynamical processes. 

Following RMT analysis of adjacency matrices of networks introduced in our pre- 
vious work [21], we investigate random matrix properties of network Laplacians. In 
this paper we only concentrate on the spacing distribution of the eigenvalues of the 
Laplacian, which explains short range correlations among the eigenvalues. Further in- 
vestigations would involve more sensitive analysis like A3 statistics to understand the 
long range correlations |36j . 
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